Discrete single-photon quantum walks with tunable decoherence 
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Quantum walks have a host of applications, ranging from quantum computing to the simulation 
of biological systems. We present an intrinsically stable, deterministic implementation of discrete 
quantum walks with single photons in space. The number of optical elements required scales linearly 
with the number of steps. We measure walks with up to 6 steps and explore the quantum-to-classical 
transition by introducing tunable decoherence. Finally, we also investigate the effect of absorbing 
boundaries and show that decoherence significantly affects the probability of absorption. 



The random walk is a fundamental model of dynami- 
cal processes that has found extensive application in sci- 
ence. The quantum walk (QW) is the extension of the 
random walk to the quantum regime [TJ [2]. Here, the 
classical walker is replaced by a quantum particle, such 
as an electron, atom or photon, and the stochastic evo- 
lution by a unitary process. A key difference is that the 
many possible paths of the quantum walker can exhibit 
interference, leading to a very different probability dis- 
tribution for finding the walker at a given location. 

An important motivation for work on QW's has been 
their application to quantum computation: not only were 
they instrumental to Feynman's original quantum com- 
puter (as the clock mechanism [3]), but it has since 
been shown that they represent a universal computa- 
tional primitive [4j [5] and have inspired novel quantum 
algorithms [6]-[8]. QW's have also been used to analyze 
energy transport in biological systems [9] [10] . 

Despite a few early experimental demonstrations [TT1 - 
fT4] , experimentalists have only recently begun to develop 
the level of control over single quantum particles required 
to implement discrete-time QW's, leading to demonstra- 
tions with neutral atoms in position space [15] . ions in 
phase space p~6j [17], and single photons in time [T8] . 
Continuous-time quantum walks have different outcomes, 
applications, and experimental challenges, see, e.g., [T9] . 

In this work, we present discrete-time QW's of sin- 
gle photons in space. Our approach is robust, due to 
the use of intrinsically stable interferometers, yet highly 
versatile — enabling control over every operation at every 
step of the walk. Of particular interest is the ability to in- 
troduce a controlled amount of decoherence, which we use 
to explore the quantum-to-classical transition. Besides 
being of fundamental interest, decoherence in QW's can 
improve the performance of certain computational proto- 
cols [20] and is particularly important in their application 
to describing energy transport [9] [10] . Finally, we inves- 
tigate the effect of introducing absorbing boundaries into 
the walk, as theoretically investigated by [2T| [22], 

The simplest random walk occurs on a one-dimensional 
lattice. The particle begins at one site and each step 



of the walk consists of a move to a neighboring site 
on the left or right, determined by the outcome of a 
coin flip. In the analogous QW, the coin is represented 
by a two-level quantum system whose orthogonal lev- 
els we will designate \H) (horizontal) and \V) (vertical). 
Each step of the QW starts with an analogue of the 
coin flip: a unitary coin operator C is applied to the 
coin space. An unbiased, or Hadamard, coin operator 
transforms the coin so that \H)^\D)=(\H) + \V))/V2, 
\V)^>\A) = (\H)-\V))/y/2. Following each coin operation 
is the shift operator, 

S = X) \j - 1)01 ® + \j + ® \V)(V\, (1) 
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which moves the particles to one of the neighboring lat- 
tice sites, conditional on the quantum coin state. There- 
fore, the operation W = SC makes up a single step of 
the QW, and a walker in an initial state is found in 
the state W N \ip} after N steps. 

The quantum and random walks can be considered 
the extremes of a spectrum, with pure quantum evolu- 
tion turning into classical evolution if there is sufficient 
decoherence [23U25] . In addition, all the intermediate 
walks are special cases of the broad category of quantum 
stochastic walks [26] . Various mechanisms of decoherence 
have been studied [25 ; here we consider pure dephasing, 
since this corresponds to the decoherence mechanism in 
our experiment. The system is described by a density 
matrix p which at each step undergoes the evolution 

PiV+i = (1 - q)Wp N Wi + q^KiWp N WiKl (2) 

i 

where the Kraus operators Ki = correspond to pure 
dephasing. The parameter q is the probability of a de- 
phasing event occuring at each step. If q = 0, the walk 
is a pure QW, while q = 1 reproduces the random walk. 
Note that to observe the quantum-to-classical transition, 
gradual decoherence must be applied at each step, not 
just to the initial coin state [23H25] , 

Our implementation of the discrete QW on a line rep- 
resents a polarization analogue of the originally pro- 
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FIG. 1: Experimental schematic, (a) Photons created via 
spontaneous parametric down- conversion in a PPKTP crys- 
tal are injected into free-space mode 0. Arbitrary initial coin 
(polarization) states are prepared by the first polarizing beam 
splitter (PBS) and waveplate combination. Six pairs of coin 
(Ci) and shift (Si) operators implement six steps of the walk. 
Coincident detection of photons at detector D2 and Dl (4.4 ns 
time window) herald a successful run of the walk, (b) Details 
of our optical mode numbering convention for the first two 
steps. The dashed lines trace out one of the interferometers 
which are used to align the quantum walk, (c) A relative an- 
gle between two beam displacers reduces the recombined pho- 
ton's temporal (At) and spatial (Ax) mode overlap, thereby 
implementing tunable decoherence. 



posed linear-optical version of the Galton board based 
on beamsplitters and phaseshifters [12j [27]. A simi- 
lar, polarization-encoded setup has been proposed for 
cube polarizing beamsplitters [28 . Figure []Ji) shows a 
schematic of the experiment. 

Pairs of single photons were created via type-II sponta- 
neous parametric down- conversion in a, nonlinear, potas- 
sium titanyl phosphate (PPKTP) crystal. This crystal 
was pumped by a 5 mW diode laser centered at 410 nm 
and emitted orthogonally polarized photon pairs with 
a wavelength of 820 nm and a FWHM bandwidth of 
0.6 nm. The pairs were separated at a polarizing beam- 
splitter; one photon from each pair served as a trigger, 
while the second photon was launched into the QW. At 
an average heralded photon rate of ^20000 s _1 , the mean 
longitudinal distance between two photons was about 
250000 times longer than the setup length of ~60 cm. 
The probability of randomly creating more than one pho- 
ton pair simultaneously was ~9 • 10 -5 , i.e., only one pho- 
ton was in the setup at any given time. 

Quantum coin states were encoded in the polarization 
\H) and \V) of the input photon. Throughout our ex- 
periment, the initial coin state was set to left-circular 
polarization, \L) = (\H)-\-i\V))/y/2, using a quarter- and 
a half- wave plate, leading to symmetric probability dis- 
tributions. For the results presented here, the remaining 
coin operators C were Hadamards, realised with half- 
wave plates set to 22.5°. We can, however, prepare arbi- 
trary pure input states as well as arbitrary coin operators 



for each step with suitable wave plate settings. 

The lattice sites of the QW were represented by lon- 
gitudinal spatial modes. The shift operator S acting on 
these modes was implemented by a 27 mm long, birefrin- 
gent calcite beam-displacer. The displacers had a clear 
aperture of 20x10 mm 2 and were mounted on manual, 
tip-tilt rotation stages with a resolution of 217 |irad/5° 
turn. The optical axis of each calcite prism was cut so 
that vertically polarized light was directly transmitted 
and horizontal light underwent a 2.7 mm lateral displace- 
ment into a neighboring mode. Lattice sites were, typical 
for discrete walks on a line, labeled so that there were odd 
sites at odd time steps and even sites at even times. 

The first two steps of the QW are shown in detail in 
Fig. [T]o) . The spatial modes after step 1 were recom- 
bined interferometrically in step 2. A series of steps then 
formed an interferometric network, Fig. [IJi). We aligned 
this network iteratively, to a single interferometer per 
step. For example, the second beam displacer was aligned 
to maximize the interference visibility of the interferom- 
eter in Fig. [TJd): the state \D) was input in mode and 
the beam displacer rotated to maximize the overlap of the 
output mode with \D). The third displacer was then 
aligned to the second and so on. We reached interference 
visibilities of typically ~ 99.8% per step. 

The photons emerging in the N+l spatial modes at 
the output of an TV-step QW were coupled into an opti- 
cal fiber and subsequently detected by a single-photon 
avalanche photodiode, in coincidence with the trigger 
photon. We measured the probability distributions se- 
quentially, translating the fiber coupler between the in- 
dividual modes using a manual translation stage. 

The measured probability distributions for detecting 
the photon at a given site, for 1 to 6 steps, are shown in 
Fig. . The experimental data are in excellent agree- 
ment with theory, with an average Li-norm distance, 
d = \ E* \Pi XP ~Pi h \i of 0.031 for the coherent and 0.017 
for the decohered walks. The quality of our data degrades 
somewhat for a higher number of steps, largely due to 
nonplanar optical surfaces, which caused small relative 
phase shifts between the multiple interferometers. The 
decohered walk is insensitive to phase errors and there- 
fore better agrees with theory. 

Our scheme has several advantages: first, the inter- 
ferometric network is inherently stable. The transver- 
sal mode-match is fulfilled because two beams emerging 
from one displacer will always be parallel, independent of 
small deviations in the optical alignment. The stability 
and alignment procedure of the QW grid are facilitated 
by the fact that the N interferometers between steps N 
and (N+l) are formed between only 2 optical compo- 
nents. Our setup — even though it is interferometric — 
does therefore not require active phase locking. Secondly, 
our system scales well, with the number of optical com- 
ponents increasing as 2 N (as opposed to (N 2 +N)/2 in 
[I2j[27]), and exhibits low optical loss of ~ 1% per step. 
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FIG. 2: Probability distributions for successive steps of the (a) quantum and (b) fully decohered (classical) walks up to the 
sixth step. Dashed lines show experimental data and solid lines show theoretical predictions. Probabilities are obtained by 
normalizing photon counts at each position to the total number of counts for the respective step. The insets show horizontal 
scans across the walk lattice for the five-step quantum walk (coupled into single-mode fiber) and decohered random walk 
(multi-mode fiber), respectively, (c) Normalized standard deviation of the probability distribution for quantum (black circles) 
and classical (red triangles) walks for 1 to 6 steps. Lines show the theoretical values; error bars are smaller than symbol size. 



The remaining obstacle to scalability are nonideal opti- 
cal components, a problem that can be alleviated with 
careful manufacturing or the use of shorter displacers. 

A unique feature of this setup is that tunable deco- 
herence can be introduced by intentional misalignment 
of QW steps, Fig. Setting a non-zero relative an- 

gle between neighboring beam displacers leads to both 
a temporal delay and a transversal mode mismatch be- 
tween interfering wave packets. Because the coincidence 
time window was much longer than the temporal shift we 
essentially integrated over the timing information, which 
corresponds to dephasing, cf. Eq. [2] Similarly, we traced 
over the spatial mode information by coupling the pho- 
tons into a multi-mode fiber — as opposed to the single- 
mode fiber used for the coherent walks shown in Fig.[2^). 
In practice, this reduced the interference contrast in all 
interferometers in the respective step. The QW was fully 
decohered [q=l in Eq. [2] when the interference visibility 
in each individual step reduced to 0, which occured at 
a relative angle of 10.5° in our experiment. Figure J2Jd) 
shows the experimental results given by our system at 
full decoherence for steps 1 to 6. The probabilities — as 
expected for a classical random walk — follow a binomial 
distribution around the origin. 

A distinguishing feature of an ideal QW is the speed 
at which the walker traverses the line. In particular, the 
standard deviation of the QW is proportional to the num- 
ber of steps and not, as for the classical walk, its square 
root [2 j. This has been exploited to design quantum- 
walk-based search algorithms that exhibit a Grover-like 
quadratic speedup [7 . The measured standard devia- 
tions for both our quantum and fully decohered walks are 
shown in Fig. [2^) . The results show very good agreement 
with theory: the fully decohered walk spreads diffusively, 



while the quantum walk spreads ballistically. 

Tunable decoherence enabled us to investigate the 
quantum-to-classical transition for a 5-step walk. By ap- 
plying Eq. [2] to a two-step walk, see Fig. [TJd), we calcu- 
lated the interference visibility in output mode after the 
second beam-displacer as a function of the decoherence 
q. We then adjusted the relative angle, Fig. [lb), between 
beam displacers to a target visibility. Figure|3] shows the 
resulting probability distributions, compared to theory, 
Eq. [2] Note the interesting nonlinear dependences of the 
probability distributions on q. 

Finally, we demonstrate another qualitative difference 
between classical and quantum walks, by incorporating 
absorbing boundaries. While a classical walker is eventu- 




FIG. 3: Transition of a 5-step quantum walk to a classical 
random walk. The decoherence parameter q for a given data 
set was obtained from a least-squares fit to theory, Eq. [2] 
All error bars are smaller than symbol size; the shaded area 
represents the theoretical probability distributions. 
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Step 

FIG. 4: Transmission probability of quantum (black circles) 
and classical (red triangles) walkers with an absorber located 
at position —1 and initial coin state \L). The transmission was 
obtained as the ratio of the number of transmitted photons 
measured with absorbers in place to the number measured 
without them. Error bars are smaller than symbol size. The 
insets show the fifth step walk with an absorber at position —1 
(red columns) compared to the original walk (blue columns) 
for the quantum (QW) and classical case (CW). 

ally absorbed in the presence of an absorbing boundary, a 
quantum walker escapes with probability 1 — 2/tt [2TJ [22] . 
A difference between the two exit probabilities first oc- 
curs after 5 steps, making it experimentally accessible 
with our current setup and providing a novel way of char- 
acterizing the degree of coherence in the walk. Absorb- 
ing boundaries were implemented using beam blocks in 
every spatial mode —1. Figure [4] shows the measured 
single-photon transmission probabilities in the quantum 
and fully decohered (classical) cases. 

The most compelling features of our scheme are the 
ability to add tunable decoherence to a QW and the 
fact that every individual lattice site is fully accessible 
at any given time step. Future work could be to im- 
plement random or position-dependent coin operators 
to study walks on random environments [29], inhomo- 
geneous walks [30] and topological insulators [31 . We 
could also add another walker on a separate line, using 
the existing setup with a vertical offset between input 
beams, or on the same line, with two or more photons 
launched in the same (or neighboring) spatial modes. 
This would allow exploration of entangled QW's [32], as 
high-quality polarization-entangled photons can be rou- 
tinely produced at high rates [33] . Finally, the setup can 
be used to prepare photon-number and path-entangled 
states across a large number of modes [34] . 
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